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Recolliding electrons are responsible for many of the interesting phenomena observed in the inter- 
action of strong laser fields with atoms and molecules. We show that in multielectron targets such 
as Ceo a new important recollision pathway opens up: the returning electron may excite collective 
modes even if the laser frequency is far off-resonant. We formulate a simple analytical theory which 
predicts that the recollision-induced excitation of collective modes should dominate over the "usual" 
harmonic generation yield at 800 nm wavelength. In this case the tomographic imaging of complex 
multielectron systems may be obscured. We employ a time-dependent density functional model of 
C60 and show that with increasing laser wavelength the dynamics becomes more and more single 
active electron-like, suggesting that long wavelengths are to be preferred for imaging purposes. 

PACS numbers: 33.20.Xx 36.40.Gk 31.15.ee 



I. INTRODUCTION 

A typical interaction scenario in strong field laser atom 
or molecule interaction involves three steps: (i) the re- 
moval of an electron from a target (ionization), (ii) mo- 
tion of this electron in the continuum, and, possibly, (hi) 
a recollision with the "parent" atom or molecule if step 
(i) occurred at a time such that the laser field drives 
the electron back. The recollision in the third step is 
responsible for the plateaus in photoelectron and high 
harmonic spectra, and nonsequential multiple ionization, 
corresponding to the three pathways (i) scattering in the 
presence of a laser field, (ii) recombination and emission 
of a photon, and (iii) laser-induced collisional ionization 
(see, e.g., [3, 0,01 for reviews). 

Structural information about the target is encoded in 
both photoelectron and harmonics spectra. Hence, be- 
sides the potential of high order harmonic generation 
(HOHG) as an efficient source of short wavelength ra- 
diation and attosecond pulses Q, the so-called "tomo- 
graphic imaging" of molecular orbitals (see [H and [j| 
for a review) has attracted considerable attention. It is 
clear that whatever is "imaged" in this procedure is sup- 
posed to be representation- independent, i.e., should not 
depend on the basis in which one expands the multielec- 
tron wavefunction. This requirement is difficult to fulfill 
within the simple and commonly adopted single active 
electron approximation (SAE) 

In this work we study the recollision dynamics and the 
emitted radiation for the case of the Ceo fullerene, which 
is an example for a multielectron system displaying col- 
lective modes and an interesting dynamics when exposed 
to fs laser pulses 0,(1] (other such systems are, e.g., metal 
clusters or biomolecules). The laser frequency is kept well 
below the surface and volume plasmon frequency of Cgo 
so that only the recolliding electron may excite the col- 
lective modes efficiently but not the laser itself. In the 
context of "orbital imaging" it is vital to know whether 
the structural information encoded in the HOHG spectra 
is "contaminated" by emission at collective frequencies. 



In other words, we are interested in the relative efficiency 
of the collective response with respect to the "standard" 
harmonic generation. 

The outline of the paper is as follows. In Sec.|TT]the Cgo 
jcllium model we use in the time-dependent density func- 
tional theory (TDDFT) calculations is reviewed, and its 
collective modes are identified. In Sec. IIIII HOHG spec- 
tra are presented for three different wavelengths, ranging 
from the typical 800 nm up to 3508 nm. The transition 
from the linear to the nonlinear excitation regime is dis- 
cussed, enhancements in the dipole spectra due to plas- 
mon excitation are evidenced, and their origin is investi- 
gated. In Sec. IIVI we compare the TDDFT results with 
the predictions of a simple, SAE Lewenstein-like model 
of HOHG from Cgo- Section [V] is devoted to an analyti- 
cal model which takes collective modes into account and 
enables us to predict the relative efficiency of harmonic 
emission due to recollision-induced plasmon excitation 
(RIPE) with respect to standard harmonic generation. 
Finally we conclude in Sec. I VII 



II. MODEL 



The Ceo fullerene is modelled using density functional 
theory (DFT) employing a jellium potential for the ionic 
background of inner and outer radius Ri, R , respectively, 
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FIG. 1: (color online). Net KS potential (black, squares), 
total density (red, diamonds), wavefunctions of the lowest 
KS orbital and the HOMO (orange, crosses and triangles, 
respectively). The a- and it- levels are indicated. Density and 
wavefunctions are scaled to fit into the plot. 



solution of the time-independent KS equation 
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FIG. 2: Dipole response of the Ceo model system. Narrow 
lines (single-particle transitions) on top of two broad struc- 
tures (surface or Mie plasmon umic — 0.7 and volume plas- 
mon Up ~ 1.4) are observed. The Mie plasmon corresponds 
to homogeneous dipole-like oscillations of the electron density 
with respect to the ions. The volume plasmon (in general a 
breathing mode) is visible in our dipole spectra since it con- 
tains a nonvanishing dipole component. The dipole strength 
is normalized such that its integral equals N = 250. 



ejtyj) = (T + V + Va + VxcMj) 



(2) 



A. Collective modes lumic and lj p 



yields the ground state configuration from which we start 
the propagation. Here, \tpj), j = 1 . . . N are the N KS 
orbitals, ej are the KS orbital energies, T is the single- 
particle kinetic energy operator p 2 /2, 



v H = I dV n(r>) 
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is the Hartree potential, 



V xc (r) 
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(4) 



is the exchange-correlation potential in exchange-only lo- 
cal density approximation (LDA), and 



(5) 



is the electron density. The N — 250 KS electrons lead to 
a spin-neutral, closed-shell ground state of spherical sym- 
metry. More precisely, we obtain 200 er-electrons (with- 
out node in the radial wavefunctions) and 50 7r-electrons 
[with one node in the radial wavefunction located close 
to the Ceo-radius R = (Ri + R )/2 = 6.7]. The free pa- 
rameter Vo = 0.68 is used to adjust the KS energy of 
the highest occupied molecular orbital (HOMO) to the 
ionization potential of Cqq, -£homo = Ip — 0.28. The 
HOMO of our model is a 7r-orbital of angular momentum 
quantum number I — A. Figure [1] illustrates and summa- 
rizes the ground state configuration from which we start 
the time-dependent calculations. 



In order to characterize the collective response of the 
model Cf;o we apply the real-time method proposed in 
Ref. To that end we solve the time-dependent KS 

(TDKS) equation [13 

[ Wt^ J {t)) = [ T + V + V ^) + V k + V ^o (*)> (6) 



with 



V I (t) = A(t)-p 



(7) 



where A(i) = Ae z Q(t) is a vector potential describing 
a 5-like electric field E(t) = —dA/dt = AS(t) in dipole 
approximation. From the Fourier-transform of the dipole 



d z (t) = / d 3 r zn(r,t) 



(8) 



the spectrum S(tu) = \d z (tu)\ 2 is calculated. Figure [5] 
shows that the linear dipole response consists of several 
narrow lines (single-particle transitions) that sit on top 
of two broad structures (the surface and volume plasmon, 
respectively). Closer inspection shows that transitions of 
the type ai — > 7r(£ ± 1), nl — * <j{1 =F 1) contribute to the 
surface (or Mie) plasmon WMie and transitions between a- 
states and (initially unoccupied) 5-states (with two radial 
nodes) to the volume plasmon w p . 



III. RESULTS 

In this Section we shall present and discuss our results 
for dipole spectra S(u>) of our model Cgo when exposed 
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to Gaussian and trapezoidal laser pulses of various peak 
intensities and wavelengths. 



A. From linear to nonlinear plasmon excitation 

We solved the TDKS equation ^ for Gaussian pulses 
with a vector potential of the form 



Alt) = - — sin(W) e~ a '® 



with 



a(t) = 3 



(9) 



(10) 



There are 0.278 n cycles within the FWHM of the Gaus- 
sian pulse (with respect to the electric field or the vec- 
tor potential) centered around t = (2n/ui)(n/2). We 
started the simulation from the ground state at t = 
and stopped at t = {2ir/uj\)n with n = 8. 

Figure [3] shows the transition from the linear to the 
nonlinear regime. At very low field amplitude (E = 
0.0025 and 0.005 at 2280 and 800 nm, respectively) the 
dipole spectra display replicas of the linear response pro- 
file on a very low level, depending on the bandwidth of 
the applied laser pulse. Upon doubling the field ampli- 
tude (E = 0.005 and 0.01 at 2280 and 800 nm, respec- 
tively) the signal in the dipole spectrum is quadrupled, 
as expected in the linear regime. The corresponding val- 
ues of intensity are 0.9 • 10 12 W/cm 2 and 3.5 • 10 12 W/cm 2 . 
However, with further increasing laser intensity, plateaus 
develop and the high harmonic-signal increases rapidly 
over a wide frequency range. One may argue that this 
increase of the harmonic signal is just due to the stan- 
dard harmonic generation mechanism while the collective 
response is still within the linear regime and thus not 
visible at higher laser intensities. The next Subsection 
is hence devoted to identify plasmon enhancements and 
their wavelength dependence. 



B. Plasmon enhancements and wavelength 
dependence 



Figure 0] shows the harmonic spectra S(u>) as calcu- 
lated from the full dipole and the outermost orbital den- 
sity only ('HOMO only') for an 8-cycle, (2,4,2) trape- 
zoidal 800- nm laser pulse, i.e., with 2-cycles up and down 
ramps and 4 cycles of constant amplitude E = 0.05 
[HI]. The difference between the two harmonic spectra 
clearly indicates that not just the valence electron con- 
tributes to the emission. Enhancements by two orders of 
magnitude around frequencies at which the system dis- 
plays collective modes are visible. The standard cut-off 
known from atomic HOHG is at 3.17f7 p + |£homo| (with 
Up = E 1 1 (4cl> 2 ) the ponderomotive energy) and indicated 
by an arrow. The real cut-off, however, is extended to 
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FIG. 3: (color online). Harmonic spectra of the Ceo model 
for Gaussian laser pulses with uj\ = 0.057 (A = 800 nm, upper 
panel) and uj\ = 0.02 (A = 2280 nm, lower panel). The values 
of E are given in the plots. The linear response profile from 
Fig. [2]is included (dotted). The field amplitude E = 0.04 
corresponds to the intensity 5.6 ■ 10 13 W/cm 2 . 



higher harmonic frequencies because recombination into 
orbitals with higher ionization potentials \ej\ > |£homo| 
takes place. Note that the latter is possible without viola- 
tion of the Pauli principle (unless KS electrons are frozen 
in the respective states). An extension of the standard 
harmonic plateau in a multielectron system — presumably 
of the same origin — has also been observed in Ref. [l4[ • 

In the following we show that with increasing laser 
wavelength the emission spectra become more and more 
SAE-like in the sense that all collective response is less 
efficient than the standard harmonic generation by the 
outermost electron at the respective frequency. In the 
SAE calculations we also start from the DFT groundstate 
but freeze the potentials Va and V xc for the propagation 
of the valence KS orbital. 

Figure [5] shows that at A = 2280 nm there are still sub- 
stantial differences between the SAE-result and the full 
TDKS calculation. First, the SAE yield is higher because 
the ionization step in the three step scenario described 
above is more efficient for a frozen potential since there is 
no polarization which counteracts the laser field. Second, 
the plasmon emission included in the full result obscures 
the oscillatory structure from which structural informa- 
tion (i.e., in our case the Ceo radius and the width of the 
spherical jellium shell) could be obtained. Only in the 
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FIG. 4: (color online). Harmonic spectra of the Ceo model 
for E = 0.05, lo\ = 0.057 (A = 800 nm), and an 8-cycle trape- 
zoidal laser pulse with 2-cycles up and down ramps The 
full spectrum and the one just from the valence KS electron 
('HOMO only') are shown. The linear dipole response from 
Fig. [2] is included (shifted vertically). The vertical arrow in- 
dicates the standard cut-off 3.17U P + £homo|- 




12 3 

Frequency (a.u.) 



FIG. 5: (color online). Emission for wi = 0.02 (A = 2280 nm) 
and E = 0.03 (other parameters as in Fig. [I). The results 
from a full TDKS calculation ('full') and a SAE-simulation 
are shown. The linear dipole response from Fig. [2] is included 
(shifted vertically) . The vertical arrow indicates the standard 
cut-off. 

(extended) cut-off region full and SAE-result agree very 
well because there are no collective modes at such high 
frequencies. 

At the even longer wavelength A = 3508 nm the full 
TDKS result agrees well with the SAE result, as is shown 
in Fig. O Also the cut-off is at the expected position, in- 
dicating that recombination into states with orbital ener- 
gies |e| > |chomo| is insignificant. A closer inspection of 
the individual response of all the KS electrons shows that 
the standard HOHG generation of the HOMO KS elec- 
trons (i.e., the two spin-degenerate ones with i = 4 and 
m = 0) clearly dominates. Hence, long wavelengths are 
advantageous for imaging schemes which are based on in- 
terference structures in the HOHG spectra predicted by 
strong field-theoretical treatments [15j in SAE approxi- 
mation. However, the efficiency of HOHG also decreases 
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FIG. 6: (color online). Emission for on = 0.013 (A = 3508 nm) 
and E — 0.02. The results from a full time-dependent KS cal- 
culation ('full') and a SAE-simulation are shown. The bold 
vertical arrow indicates the standard cut-off, the three thin 
vertical arrows local minima in the envelope of the full spec- 
trum. 

with increasing laser wavelength [Hj( . The fact that the 
efficiency of the collective response decreases even faster 
is one of the main results of this work. 



C. Identifying the mechanism 

In Sec. MI Al we showed that at low field strengths 
the collective response increases linearly with the field 
strength (i.e., the signal in the dipole spectra quadrati- 
cally) while at higher intensities the standard high har- 
monic plateau develops, which is an entirely nonlinear 
phenomenon. In Sec. IIIIBl we showed that plasmon en- 
hancements are present although they decrease relative 
to the standard high harmonic plateau with increasing 
wavelength. This means that there must be some non- 
linear effect at work which is able to generate a collective 
response of comparable strength as the standard high 
harmonics. The latter are due to returning electrons 
which rccombine. The obvious guess is to attribute the 
collective response also to the returning electrons so that 
the similar efficiency of harmonic emission via the SAE 
and via the collective mechanism can be understood if 
recombination with emission of a photon and with exci- 
tation of a plasmon (followed by emission of a photon) 
are similarly efficient. In this Subsection we support the 
viewpoint that the recolliding electrons indeed excite col- 
lective modes by analyzing our numerical results in more 
detail. 

In our TDDFT simulations we use a spherically sym- 
metric imaginary potential W(r) = — iVt / o(r/i? gr id) 16 
with Wq = 100 and i? gr id the radius of the numerical 
grid. The imaginary potential serves as an absorber of 
probability density approaching the boundary of the nu- 
merical grid [k| • Usually the grid is chosen big enough so 
that only the probability density corresponding to never- 
returning electrons is absorbed and thus the imaginary 
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potential does not affect the relevant dynamics taking 
place in the interior of the numerical grid where W(r) is 
negligible. However, in order to test whether recolliding 
electrons are responsible for both the standard harmonic 
generation and the plasmon enhancements, we may ab- 
sorb probability density representing electrons of a cer- 
tain excursion amplitude z by moving the imaginary po- 
tential closer to the C^o- If the plasmon enhancements 
are due to recolliding electrons we then expect the har- 
monic signal and the plasmon signal to drop. If, instead, 
the plasmon enhancements are due to some other yet un- 
known nonlinear effect which does not require returning 
electrons, then the harmonic signal should drop while the 
plasmon signal sustains. 

Figure [7] shows dipole spectra for A = 2280 nm and 
E = 0.01 (i.e., the second highest intensity shown in the 
lower panel of Fig. ^ for two grid sizes. The excursion 
amplitude of a free electron in this case is z = E/u>f = 25. 
Hence we expect the i? gr id = 100-grid to comprise all the 
relevant electron dynamics whereas on the i? gr id = 40- 
grid some electrons will be already inhibited from re- 
turning to the Ceo because the corresponding probability 
density is absorbed. In fact, Fig. [7] shows that parts of 
the plateau are removed in the spectrum for the smaller 
grid. Only the single particle transition lines close to 
lu = 0.6 are unaffected by the absorbing boundary, show- 
ing that these transitions are not excited by recolliding 
electrons but — presumably — by multiphoton resonances. 
However, besides these resonant transitions the whole 
plateau is suppressed. We thus conclude that the return- 
ing electrons are essential for the excitation of the col- 
lective modes. This conclusion is further supported by a 
time- frequency analysis of the dipole d z (t) . To that end a 
spectral window is applied to d z (uj) . The result is trans- 
formed back, which corresponds to the spectral filter- 
ing of certain harmonics for the generation of attosecond 
pulses in experiments 0. The result is shown in Fig. [5] 
The emission follows overall nicely the classical "sim- 
ple man's theory" : the classical return-times of electrons 
with return-energy £ rct (which contribute to the emission 
of harmonic radiation at a frequency w = £ rc t + |chomo|) 
are indicated by white trajectories in the frequency-time 
plane. It is seen that the plasmon emission is correlated 
with the return of electrons. Whenever there are recollid- 
ing electrons having the right energy to excite a plasmon, 
enhanced emission is observed. Due to the large width of 
the collective resonances the emission decays before the 
next returning electron collides. 



IV. SAE LEWENSTEIN MODEL VS 
LONG- WAVELENGTH TDDFT- RESULT 

We now show that the structure in the HOHG spec- 
trum of Fig. [6] is indeed similar to what one expects from 
the strong field approximation applied to HOHG, i.e., the 
so-called Lewenstein-model [17[ ■ Within the Lewenstein- 
model the dipole expectation value for an infinite, lin- 
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FIG. 7: (color online). Same as in Fig- El lower panel, for E = 
0.01 but two different grid sizes (indicated in the plot). The 
linear response profile from Fig. [5] is included again (dotted). 




FIG. 8: (color). Logarithmically scaled contour plot of the 
time-frequency analyzed dipole emission log 10 |d z (£,a;)| 2 for 
the parameters of Fig. [5] The white lines indicate the classi- 
cal solutions of returning electrons (see text). The positions 
of the Mie surface plasmon and the volume plasmon are indi- 
cated. 



early polarized laser pulse 

E{t) = Ee z cosoj l t, E(t) = -d t A(t) (11) 
is given by 

/•co /9_\ 3 / 2 

d^(t) = i ^ dr^-J v*Mt,T)+A(t)} (12) 

x exp[— iS(t,r)]E cos[uJi(t — r)] 
X(i z [k(t,T) + A{t- t)} +c.c. 

where r is the travel-time of the electron between ioniza- 
tion and recombination, 

»z(pz) = (p*N*o>, (13) 
k(t, t) is the saddle-point momentum 

- COSLU]t — COS Wl ft — t) 

k(t,r) = -E ! '-, (14) 
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FIG. 9: (color online). Harmonic spectra S(u), calculated 
from Eq. (fT2)l for I = 4, I = 0, and I — 4 but twice the radius 
R. The laser parameters are the same as in Fig. [()] [wi = 0.013 
(A = 3508 nm) and E = 0.02]. R = 6.7 and A = 1.4 was 
used. The three thin vertical arrows indicate local minima in 
the envelope of the spectrum for I = 4. 



and S(t, t) the saddle-point action 



S(t,r) = (C/p-£HOMo)l" 



2{7 P (1 — cos u>\t) 



U p C(t) cos[(2t - t)u{\ 



(15) 



with 



and zero otherwise. Here, A is half the thickness of the 
Cgo-shell, i.e., A = (R — Ri)/2 = 1.4. Assuming further 
\p z A\ <C 1 and, e.g., £ = we obtain 

Hz{Pz) ~ \{smp z R-p z Rcosp z R+p 2 z R 2 smp z R) (17) 
PS 

and a similar but more lengthy expressions for £ = A. 
One clearly sees that structural information (i.e., the Ceo- 
radius R) is 'encoded' in /j, z (p z ). If the approximation 
|p z A| <C 1 is not made, also information about the shell 
thickness 2 A is included in the matrix element /j, z (p z ). 

Figure [9] shows the harmonic spectra obtained from 
the Fourier-transform dz (uj) of Eq. (|12[) for £ = 4 and 
£ = and the laser parameters of Fig. [5J The positions 
of the minima in the envelope of the HOHG spectra de- 
pend on the initial £ quantum number and the Ceo radius 
R. In order to illustrate this dependency the spectra for 
£ = 0, £ = 4, and £ — 4 but with the radius doubled 
are shown. The minima indicated in the £ = 4-spectrum 
by vertical arrows may be compared with those of the 
time-dependent DFT result in Fig. [51 The latter are at 
uj ~ 0.5, 0.95, and 1.7. The arrows in Fig. [5] are at 
uj ~ 0.52, 1.05, and 1.9, which is in reasonable agree- 
ment. Note that the agreement would be worse if one 
attempted to compare with the £ = 0-spectrum, let alone 
with the spectrum for I = 4 and doubled radius, which 
is qualitatively different since there is at least one more 
pronounced minimum in the envelope. 



C(r) = sinter sin 2 (w!r/2). (16) 

U)\T 

For a derivation of (TT2")) the reader is referred to the orig- 
inal work in Ref. pTj . 

The target-dependence of the HOHG spectra enters 
in (fT2")) via the initial state ^> through the ionization 
and recombination matrix elements u z \k(t, r) + A(t — r)] 
and ^[k(t, t) + A(t)], respectively [13]. We assume an 
initial state of the form ^>o(r) = $o(r)Y£o(9,ip)/r with 
Ytm(0, a spherical harmonic and model the valence ir- 
orbital using a radial wavefunction $o(r)/r = (2A) -1 / 2 
for R - A < r < R, -(2A)" 1 / 2 for R < r < R + A, 



V. ANALYTICAL MODEL FOR HARMONIC 
GENERATION INCLUDING 
RECOLLISION-INDUCED COLLECTIVE 
EXCITATIONS 

For systems with a single active electron, harmonic 
spectra are usually analyzed using the strong field ap- 
proximation (or Lewenstein model) [U Gil, as we did 
in the previous section. Emission into a mode with fre- 
quency u> and polarization e\, A = 1, 2, by a system with 
only a single active electron can also be described by the 
amplitude 



//>oo />t 
d 3 p / di<*o,Wl^|0 Ul A,p) / dt'(p\V\* ). (18) 
J — oc J — oo 

Here, |^o) is the single-electron ground state, \p) is the Volkov state of drift momentum p, \n^,\) is the Fock state 
of the harmonic radiation field with n photons in the respective mode (n — 0, 1 in our case) and V3 d and V are the 
interaction operators coupling to the radiation and the laser field, respectively, 



= -i\l 2 -^r ■ e A a+ x a"* (19) 



V = E(t)z, 



(20) 
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with a + the photon creation operator and V the quantization volume. The amplitude (fl8|) describes an electron 
which is lifted from the ground state to a Volkov state by the laser field at time t' and emits a harmonic photon 
upon recombination at time t > t' . The harmonic spectrum is given by the square modulus of (fl8|) and appears to 
be virtually identical to the spectrum found from the dipolc (|12p ,T9j . In the dipole approximation we use here the 
wavevector K of the emitted photon does not appear in the amplitude l|18p. 



Now we introduce a similar amplitude which accounts for the collective modes: in addition to the pathway described 
by (fT5)) the recombining electron may excite collective modes which then relax upon emission of a harmonic photon. 
The amplitude for such a process reads 



/OO i> i>t i>t 

d^l^AlF+d |0„,A,Li) d 3 p / dt'i^LjplO^p) dt"(p|V|# >. (21) 
. . -OG <J J — OO J — OO 



Here, and \Lj) are the ground and the L-th excited 
state of a collective mode, labelled by j (e.g., surface or 
volume oscillations) . The interaction energy between the 
electron and the residual electron cloud is described by 
the operator U. The amplitude l|2"Tj) is a straightforward 
generalization of the Lewenstein model to the case when 
collective modes can be involved in the emission process. 

To evaluate the collective amplitude (j2Tj) a certain 
model for the description of the collective modes and 
their interaction U with the active electron is required. 
In order to estimate the relative contribution of the SAE 
and the collective pathways to the radiation spectrum we 
use a simple model which takes collective degrees of free- 
dom into account as two noninteracting harmonic oscilla- 



tors with eigenfrcquencies u>Mie — 0.7 and 



1.4 (i.e., 



the surface and the volume plasmon in the Cgo-niodel 
above). The respective widths of the plasmons are taken 
as r M i e ( p ) w 0.2. The main physical mechanism which 
generates these widths is a coupling between collective 
and single-electron degrees of freedom. This can also 
be interpreted as collisionlcss or Landau damping of col- 
lective modes in a finite system [2(| [2l|. To obtain an 
explicit form for the interaction operator U we employ a 
rigid sphere model (RSM) in which the electron cloud is 
treated as an incompressible homogeneous sphere which 
may oscillate around its equilibrium position. Note that 
on the level of modelling in this Section it does not matter 
whether we consider a homogeneous sphere or a spherical 
shell. Within the RSM the interaction operator has the 
form 



U(r,X) 



(N-l) 
R 



(r-xy 



R 



2R 2 



X\ 



\r-X\<R, 
\r-X\ > R, 



(22) 

where N — 240 is the number of electrons, r is the active 
electron's position and X is the center-of-mass displace- 
ment of the electron cloud. Because of the relatively high 
energies of the plasmons only the first excited collective 
states are relevant in the sum over L in (12111. For a first 



excited state X ~ 1/ y/(N — l)cjMic(p) <C i? so that with 
high accuracy (e.g., taking N = 240 and wmic = 0.7 one 
estimates X ~ 0.08) one may simplify (f2"2"|) keeping only 
the linear term with respect to the center-of-mass dis- 
placement X: 

U(r,X)^U (r) + ^r. X {^ /r3j ^ . (23) 

Next, we assume that the electron excursion amplitude 
in the laser field z — E/uf is less than or comparable to 
the cluster size R. Then, with reasonable accuracy, we 
may use 



N-l 
U(r,X)^^ r r-X 



(24) 



instead of (|23| . Within this approximation an explicit 
relation between the amplitudes (|2~lj) and (Tj"5|) can be 
derived. To this end we first evaluate the emission matrix 
element in (fT8|) . 



2irui 



<*o, WlKtdl<W,P> =- i V— (*o|r-e A |p) c iwt . 

(25) 

A similar procedure for the emission matrix element in 
(12"T|) yields 



(0Mie(p)) lu^A l^ad 1 l-Mic(p)) 



(26) 



' 7r(AT- 



e z ■ e A c 



— WMie(p) )t — r Mie ( p ) t/2 



Here we use d the f act that for the harmonic oscillator 
(0\z\l) = 1/V2MQ with M = (N-l)m and fl = uJ Mic{p) 
in our case. Also we take into account that the oscillator 
is exited along the polarization direction given by the 
unit vector e z . 

Rearranging the time-integrations, the amplitude (|2ip 
can be also written as 
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M con (u>, A) = / d W dt ' ^o, lMic(p)l^|0 M ic(p),p) / dt (0 M ic( P ), lu,\\V+ d K,\, l M ic(p)) / dt" (p\V\* ). 

Mil J J — oo J t' J — oo 

(27) 



j Mie,p 

Now the inner integral over t can be evaluated explicitly using (|26| . The result reads 



dt (0 Mi e( p ), lw.xI'^.^lOa, a, l M io( P )) = -iW ~r; ^ — - — 7^ r, r- (28) 

V VwMie(p) r M io(p)/2 - l(w ~ ^Mie(p)) 



Finally, using the standard expression for the coordinate matrix element of the harmonic oscillator (0|z|l) = l/\/2MO 
and (1241) one obtains for the first matrix clement in (l27l) 



(*o,lMic(p)l^|0 M ie(p),p) = 1 (*oNp) (29) 

" y ^Mie(p) 

Collecting Eqs. ([25]) — (|29|1 . we may express the amplitude (f2Tj) via (fT8|) as 

Mcou = i^— i ( ^ c r + r 1 M SAE . (30) 

2R A I r Mie /2 - i(w - WMieJ r p /2 - 1(U - LJ p ) 



Equation (|5D|) shows that collective modes may lead to 
enhancements in the HOHG spectrum around the respec- 
tive plasmon frequencies. For the plasmon enhancements 
to be detectable |M coU | 2 > |M SAE | 2 should hold. For the 
ratio of collective to SAE HOHG efficiency we obtain 



|M coll | 2 



|M SA e(^ = w Mio(p) )| ; 



N — 1 



R 3 uj m 



ie (p) Mie(p) 



(31) 



For N = 240, R = 6.7, r M ic = T p = 0.2 the ratio (ST]) is 
above 10 for the surface and about unity for the volume 
plasmon. 

The ratio (f3"Tj) does not depend on the laser parame- 
ters anymore whereas in our TDDFT results we observe a 
wavelength-dependent relative efficiency of the plasmon 
enhancements. With increasing laser intensity or wave- 
length the electron's excursion amplitude is increasing 
and the approximation (|24p for the interaction between 
the active electron and the electron cloud becomes in- 
valid. Without the assumption of small excursion ampli- 
tudes (as compared to the cluster radius) a simple rela- 
tionship of the type (|5D|) cannot be established. Quali- 
tatively it is quite obvious, however, that with increas- 
ing excursion amplitude distances r ~ R [for which (|23p 
is sizeable] contribute less and less to the spatial ma- 
trix element (j2"!)|) . As a consequence the standard single- 
electron HOHG spectrum dominates for E/cuf ^> R. In 
fact, E/lu? = 15.4, 75.0, and 118.3 in Figs, gj and® 
respectively, supporting our statement. 

The results ([5TJ| and (|3"Tj) were derived making several 
approximations besides the one of small excursion ampli- 
tudes. For example, the surface and the volume plasmons 
were treated as independent. This makes sense if they 
are well separated from each other, i.e., |a->Mie - w p | > 



(rMic + Tp)/2, which is actually not fulfilled in the case 
of Ceo- Another simplification was that we applied the 
RSM for the description of the electron cloud. Within 
this model the volume plasmon simply does not exist. In 
a more realistic description one should use two different 
interaction potentials instead of (|2"2"|) alone, which will 
lead to two different coefficients in (|3~T1). 



VI. CONCLUSIONS 

In conclusion, we predict a new recollision effect in 
the interaction of strong laser fields with multi-electron 
systems. Besides the usual high-order harmonic gen- 
eration the rccolliding electron may excite collective 
modes instead of emitting its energy directly as a har- 
monic photon. Via the recollision mechanism collective 
modes can be excited even if the incident laser is far 
off-resonant with the plasmon frequencies. Using time- 
dependent density functional theory we have studied the 
wavelength-dependence of the process in the case of Ceo- 
With increasing laser wavelength the dynamics becomes 
more and more single active electron-like. Experiments 
employing imaging techniques based on recolliding elec- 
trons are hence more likely to reveal clean structural in- 
formation if sufficiently long wavelengths are used. 
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